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I. INTRODUCTION 

In the standard model [1], the mixing of quarks origi- 
nates from their interactions with the Higgs field causing 
the quark mass eigcnstatcs to be different from the quark 
flavor eigenstates. Constraints on the mixing phases are 
obtained from measurements of the decays of neutral 
mesons. Decays to final states that are common to both 
partners of a neutral-meson doublet are of particular im- 
portance. The interference between the amplitude of the 
direct decay and the amplitude of the decay following 
the particle-antiparticle oscillation may lead to a CP- 
violating asymmetry between decays of mesons and an- 
timesons. The decays — > J/tpX, where X stands for 
a pair of charged kaons or pions, are a sensitive probe for 
new phenomena because the CP- violating phase that ap- 
pears in such decays is predicted in the standard model 
to be close to zero with high precision [2]. 

The first observation of the decay sequence — > 
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,7/^/2(1525), /a(1525) K+K-, was recently reported 
by LHCb [3]. The spin J = 2 assignment for the K + K~ 
pair was based on excluding pure J = 0. In this article, 
we confirm the presence of the decay B® — >• J/ipK + K~ 
with K + K~ invariant masses M(K + K~) close to 1.52 
GeV and we determine that the spin of the K + K~ res- 
onance is consistent with J = 2 and is preferred over 
,7 = or J = 1 assignments. We identify the resonance 
as the ^(1525) meson and measure the branching frac- 
tion relative to the well established decay B® — > J/i(j(j). 



II. DETECTOR 

The DO detector consists of a central tracking system, 
a calorimetry system and muon detectors, as detailed in 
Refs. [4-6]. The central tracking system comprises a sili- 
con microstrip tracker (SMT) and a central fiber tracker, 
both located inside a 1.9 T superconducting solcnoidal 
magnet. The tracking system is designed to optimize 
tracking and vertexing for pseudorapidities \r)\ < 3, where 
r] = — ln[tan(0/2)] and 9 is the polar angle with respect 
to the proton beam direction. 

The SMT can reconstruct the pp interaction vertex 
(PV) for interactions with at least three tracks with a 
precision of 40 fjjn in the plane transverse to the beam di- 
rection and determine the impact parameter of any track 
relative to the PV with a precision between 20 and 50 fj,m, 
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depending on the number of hits in the SMT. 

The muon detector surrounds the calorimeter. It con- 
sists of a central muon system covering the pscudorapid- 
ity region \r]\ < 1 and a forward muon system covering 
the pseudorapidity region 1 < \rj\ < 2. Both systems 
consist of a layer of drift tubes and scintillators inside 
1.8 T toroidal magnets and two similar layers outside 
the toroids. 



mass [1], and p is the particle momentum. To increase 
£?|? signal purity and reject prompt background, we re- 
quire the proper decay length to be greater than 200 ^m. 
The distribution of the uncertainty on the proper decay 
length peaks around 25 (im and has a long tail extending 
to several hundred microns. To remove poorly recon- 
structed events in the tail, we require the uncertainty of 
the proper decay length to be less than 100 fim. 



III. EVENT RECONSTRUCTION AND 
CANDIDATE SELECTION 

The analysis presented here is based on a data sample 
corresponding to an integrated luminosity of 10.4 fb _1 ac- 
cumulated between February 2002 and September 2011 
at the Fermilab Tevatron collider. Events are collected 
with single-muon and dimuon triggers. Some triggers re- 
quire the presence of tracks with a large impact param- 
eter with respect to the PV. The events selected exclu- 
sively by these triggers are removed from our sample. 

Candidate — > J/i/jK + K~ , J/ip — > fi + ^~ events are 
required to include two oppositely charged muons accom- 
panied by two oppositely charged tracks. Both muons 
are required to be detected in the muon chambers inside 
the toroidal magnet, and at least one of the muons is 
required to be also detected outside the toroid. Each of 
the four tracks is required to have at least one SMT hit. 
In addition, the kaon candidates are required to have at 
least two hits in the SMT, at least two hits in the central 
fiber tracker, and a total of at least eight hits in both 
detectors. 

To form B® candidates, muon pairs in the invariant 
mass range 2.9 < M(/x + ^~) < 3.3 GeV, consistent with 
J/ip decay, are combined with pairs of oppositely charged 
particles (assigned the kaon mass), consistent with pro- 
duction at a common vertex. A kinematic fit under the 
_B" decay hypothesis constrains M(fi + fi~) to the Particle 
Data Group [1] value of the J/ ip mass and the four tracks 
to a common vertex. The trajectories of the four B® de- 
cay products are adjusted according to this kinematic fit. 
In events where multiple candidates satisfy these require- 
ments, we select the candidate with the best fit proba- 
bility. The x 2 °f the fit is required to be less than 15, for 
a total number of degrees of freedom of 4. We require 
that the transverse momenta of the B® and K ± mesons 
are larger than 8 and 0.7 GeV, respectively. To suppress 
background from the decay B° — > J /ipK* (892), we re- 
quire the kaon pair to have an invariant mass greater 
than 1 GeV under the K ± tt t hypothesis. 

To reconstruct the PV, we select tracks that do not 
originate from the candidate B® decay and apply a con- 
straint to the average beam-spot position [4] in the trans- 
verse plane. We define the signed decay length of a £>" 
meson, L^ y , as the vector pointing from the PV to the 
decay vertex, projected on the B® transverse momentum 
Pt- The proper decay time of a B® candidate is given by 
t = AIb,L^ -pj (p\) where Mb s is the world-average B® 



IV. MONTE CARLO SAMPLES 

Some aspects of this analysis require information that 
cannot be derived from data. We rely on Monte Carlo 
(MC) simulated samples to derive templates of the dis- 
tributions of the signal B® — > J /ipf '2(1525) and of the 
main background components, and to derive the de- 
tector acceptance as a function of decay angles, and 
the relative acceptance for the decays B® — >• J/ip'fi an d 
B° s -> J/V/^1525). 

We use pythia [7] to generate B® mesons and 
evtgen [8] to simulate their decay. In all simulated sam- 
ples the final states are assumed to have no polarization. 
The samples have been processed by the detector simula- 
tion and the standard event reconstruction. To take into 
account the effects of the instantaneous luminosity, the 
MC samples are overlaid with data events collected dur- 
ing random beam crossings. We have generated events 
containing the decays B° -> J /ij> ^(1525), B° -> J/ip(p, 
B° -)• J/ipKZ (1430), and B° -> J/ipK$ (1430). 



V. B° -> J/tPK+R- SIGNAL EXTRACTION 

The distribution of the B® candidate mass 
M(J/if)K + K-) as a function of M(K+K~) for ac- 
cepted events is shown in Fig. 1. The data show a 
structure consistent with the decay B® — ► J / tp f 2(1525) . 
Another possibility for the observed structure near 
M{K+K-) = 1.5 GeV is the decay B° s -> J/i/>/o(1500). 
The two states have very different branching frac- 
tions to 7r + 7r~ and K + K~. The branching frac- 
tions [1] are B(/o(1500)^7T7r)=(34.9 ± 2.3)%, 
B(f (1500)^KK)=(8.6 ± 1.0)%, S(/^(1525)-> 
7T7r)=(0.82 ± 0.15)%, and B(/£(1525)-> KK=(88.7 ± 
2.2)%. If the M(K+K-) peak is due to the / (1500), 
a larger peak is expected near M(tt + it~) = 1.5 GeV. 
If the M(K+K~-) peak is due to the /a(1525) meson, 
a negligibly small peak is expected in the M(ir + Tr~) 
distribution. From the lack of significant B° s signal in 
the 7T + 7r - channel, presented in Fig. 2, we conclude that 
the observed state is not the /o(1500) meson. 

However, as observed by others [3, 9], a similar distri- 
bution may also result from decays of B° mesons where 
the J/tp meson is accompanied by a K* resonance, and 
a pion from the decay K* — > i^ ± 7r =F is assigned the kaon 
mass. 



>1.8h 

a) 
O 

C71.61 



S1.4[- 



D0 Run II, 10.4 fb 1 



1.2 



_L 



_L 



5 2 5 4 5 6 

M(^VK + K") (GeV) 



> 



o 

CM 



(A 

c 

CD 
> 
LU 




5.4 5.6 

M(^VK + K") (QeV) 




x1(T 



5.8 
(C) 



> 



O 

m 

(A 

■4— » 
c 

CD 

> 

LU 




1.4 1.6 1.8 

M(K + K") (GeV) 



5.4 5.6 

M((^ + lj,"K + K") (GeV) 



x10 J 



1.4 1.6 1.8 

M(K + K") (GeV) 



5.8 
(f) 




FIG. 1: (a) Invariant mass of the B g candidates as a function of M(K + K~) as well as (b) and (c) one-dimensional projections. 
The observed structure is consistent with a decay — > J /ipK + K~ proceeding through /2(1525) or /o(1500) mesons, (d) 
Invariant mass of the B® mesons as a function of M(K + K~) for the simulated decay B® — > J/^/2(1525) as well as (e) and (f) 
one-dimensional projections. 



There are two Kj (1430) states, degenerate in mass but 
differing in width, T = 0.109 ± 0.005 GeV for J = 2 and 
r = 0.27 ± 0.08 GeV for J = [1]. Simulated distri- 
butions of B® mass versus M(K + K~) for misidentificd 
decays B° -> J/ipK$ (1430) and B° J/^iq(U30) are 



shown in Fig. 3. In the case of J = 0, we use the full 
1=1/2 Kir elastic scattering S- wave amplitude [9] com- 
posed of -KTq(1430) and a nonresonant term. 

As seen in Fig. 3, the decays B° J/ipK}(U30) 
with J = 0, 2 can mimic the decay B® — > J /ijjf! 2 (lh2h) 1 
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FIG. 2: Distribution of the invariant mass M(tv + tt ). The solid and dashed vertical lines represent the world-average mass 
and the natural width of the /o(1500) meson [1], respectively. 



with the apparent B® mass peak position increasing with 
M(A + A~). We take this peaking background into ac- 
count by constructing templates of the mass distri- 
bution in steps of the K + K~ invariant mass of 50 MeV. 
Examples of the templates are shown in Figs. 4 and 5. 
Figure 6 shows the simulated B° — > J/-0/2(1525) signal 
in the same M(K + A - ) ranges, with fits to a sum of two 
Gaussian functions. 

Other possible sources of peaking background include 
B-meson decays to J/ip mesons accompanied by the 
J = 1 and J = 3 resonances A'*(1410) and A*(1780). 
The former decays predominantly to A'* (892). When the 
kaon mass is used for both tracks, the latter resonance 
would peak at M(A+A") greater than 1.78 GeV. 

The candidate mass distribution in the range 1.45 < 
M(K + K'-) < 1.60 GeV and | cos VI < 0.8 is shown in 
Fig. 7, where ip is the helicity angle defined in Sec. VI. 
A fit of templates for the B® — > J/ipK + K~ decay, the 
B° — » J/V>A|(1430) and S- wave A7r background and a 
linear combinatorial background, yields 534 ± 101 B® — >• 
J /ipK + K~ events. The relative rate of the S and T> Kit 
wave is constrained to the ratio of 3:2 reported in Ref. [9]. 

To extract the B® signal yield, we use the B° mass 
distribution for M(K + K~) between 1.35 and 2 GeV. We 
fit the simulated signal templates to the data, with the 
mass parameter of the core Gaussian function in each 
M{K+K~) bin scaled by a factor of 0.9982 ± 0.0008 ob- 
tained by matching the simulation and data as shown 
in Fig. 7. Figure 8 shows the B® mass fits using the 
templates for the signal B° — > J / i(j f 2(1525) and for the 
B° -» J/VA 2 *(1430) and B° ->■ J/VA^(1430) reflections. 
In these fits we allow the relative rates of the iS and 
T> Kn contributions to vary. Note that a nonrcsonant 
K + K~ component is implicitly included in the signal 
part of these fits. In addition to the peaking background 
there is a background due to random combinations and 
partially reconstructed B-meson decays, described by a 
linear function. We allow the relative normalization of 
the two Kj states to vary in each fit. The normalization 
parameters of the B® signal and background components 



are not constrained to be positive in order to obtain un- 
biased results for rates close to zero. We have conducted 
toy MC ensemble tests and we confirm that there are no 
biases on signal yield introduced by the described fitting 
procedure. 

As seen from Fig. 8, the fitted yields for the B° — > 
J/VA|(1430) and B° J/V>A *(1430) decays exceed the 
B® signal in most of the 11 subsamples, although the data 
do not provide a significant constraint on their relative 
strength. For an independent study of this background, 
we select events in the range 5.4 < M(J/i/jK + K~) < 
5.6 GeV, where the B° ->■ J/V>A}(1430) decays domi- 
nate over the B a s signal. In Fig. 9 we show the M(A ± 7T =F ) 
distribution for these events. A fit of a relativistic Brcit- 
Wigncr J = 2 resonance at a fixed mass of 1.43 GeV 
and with a floating width, over a background described 
by a second-order polynomial function, yields 3386 ± 390 
A* (1430) resonance events. This is in agreement with the 
total number of events ascribed to the A* (1430) reflec- 
tion in this mass range. The best fit result for the width 
is T = 0.162 ± 0.019 GeV that is in between the widths 
of the J = 2 and J = states. This study shows that we 
cannot establish the dominance of one background com- 
ponent over the other with the present data. Figure 10 
shows the B® signal yield versus M(K + K~ ) from data 
after taking into account the B° — > J /xpK^il^iO) and 
B° — > J/i/A'g (1430) templates and a linear background. 



VI. SPIN OF K+K~ STATE 

In Section V, we have presented evidence for the decay 
B° s ->■ J/iPK+K- in the range M(K + K~) > 1.35 GeV. 
The M(A+A") distribution peaks near 1.5 GeV. In this 
section we study the spin of the K + K~ system in the 
range 1.45 < M(K + K~) < 1.60 GeV. 

A A + A~ system can be in any natural parity state, 
J p = + ,l~,2+, etc. We consider the values J = 0, 1, 
and 2 for the spin of the observed structure. 
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FIG. 3: (a) Invariant mass of the candidates versus M(K+K~) for simulated decay B° -> J/ipK%( 1430) ; Kl -t i^K*, 
with the pion assigned the kaon mass, as well as (b) and (c) one-dimensional projections, (d) Kn 5-wave contribution [9] to 
the decay B° — > J/ipK ± ir^ 1 , with the pion assigned the kaon mass as well as (e) and (f) one-dimensional projections. 
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FIG. 6: Invariant mass of B® mesons from the simulated decay B® — > J/t/)/2(1525) for a sampling of different ranges 
of M(K + K~). The distributions are fitted with a sum of two Gaussian functions with free masses, widths, and relative 
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FIG. 7: Invariant mass of the B® candidates with 1.45 < M(K + K~) < 1.6 GeV. Curves show (a) the result of the fit allowing 
for a free B® signal yield (dashed-dotted lines) and a background composed of a 3:2 mixture of B° — > J/tpK ± TT^ 1 decays with 
the K tt^ system in the J = and J = 2 states (dashed lines) and a combinatorial component described by a linear function 
(dotted lines), and (b) the result of the fit assuming the same background model but setting the B^ signal rate to zero. 



The final state can be described by three independent 
angles. We define them as follows: Oh is the angle be- 
tween the direction of ^ + and B® direction in the J/ip 
rest frame, xp is the angle of the K + meson with respect to 
the B® direction in the K + K~ rest frame, and 4>h is the 
angle between the two decay planes, as shown in Fig. 11. 
The angular distribution for the decay of a spinless me- 
son into the spin-one meson J/ip and a meson of un- 
known spin J can be expressed in terms of H± = cos 9h, 
H-2 = cos-0, and <pn as follows [10]: 

dT ,9 

— ex |EA Jm rn#i,^r)>7 m (-#2,o)| 2 £Kfi), (i) 

where YJ™ are spherical harmonics, Aj m are complex am- 
plitudes corresponding to spin J and helicity m, and Q 
is either Hi, <pn or H 2 , and the sum extends over equal 
hclicitics of the daughter particles, m = or ro = ±1. 
The factor D is the acceptance of the event selection. Its 
dependence on the three angular variables is shown in 
Fig. 12. 

Due to limited statistics and a large background, we 
focus on the cos ip distribution obtained by integrating 
the angular distribution over cos 9jj and 4>h , taking into 
account the variation of the acceptance as a function 
of cos Oh- We extract the B® signal rate as a func- 
tion of I cos ip I by fitting the candidate mass in five re- 
gions of I cos V 1 1- The data and fit results are shown in 
Fig. 13. The resulting distribution, corrected for accep- 
tance, is shown in Fig. 14. Systematic uncertainties due 
to the shape of combinatorial background, signal model, 
and acceptance, are added in quadrature. In the region 
I cos'0| > 0.8, the large background prevents obtaining a 
reliable fit. 



For J = 0, the expected distribution is isotropic. For 
J = 1, the cosip distribution without the acceptance fac- 
tor is given by: 

oc Fi (2 cos 2 7/;) + (1-F 10 ) sin 2 ijj, (2) 

rfcos^ 

where Fiq is the ratio of the rate J = l,m = to the 
total J = 1 rate. For a superposition of J = and 2, 
with a free relative normalization, the cos ip distribution 
is obtained from 



dT 

dcosip 

v/F 20 (l-Fo)(5/4)(3 cos 2 ^ - 1) + cxp(i<S ) v^b| 
+ y (1 - F 20 )(l - F ) sin 2 ^(1 - sin 2 ^), (3) 

where F20 is the ratio of the rate J = 2, m = to the 
total J = 2 rate, and Fq is the J = fraction with 
relative phase angle So- 

Figure 14 shows that the data favor J = 2, hence the 
peak is identified with the /2(1525) meson. The fit prob- 
abilities for pure J = and pure J = 1 are 2.8 x 10~ 2 
and 9.8 x 10~ 3 , respectively. For J = 2 the fit probabil- 
ity is 0.27. The data are also consistent with a coherent 
superposition of J = and J = 2 states. With F 2 o = 1, 
we obtain the 5-wave fraction of Fq = 0.06 ± 0.16 and a 
fit probability of 0.37. 
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FIG. 8: Invariant mass of candidates for the decay B® — > J/ipK + K~ in a sampling of different M(K + K~) ranges. Each fit uses a 
template derived from the B° signal simulation, a combination of the templates for the decays B° — > J/f/'-Kj(1430), 7tj(1430) — >■ 
K ± ty^ 1 , J = 1, 2, as shown in Figs. 4 and 5, and a linear function describing the combinatorial background. The fit is performed 
in the range 5.25 < M(J/i{jK + K~) < 5.7 GeV, used to avoid the steeply falling background from multibody decays of B 
mesons at lower masses and the steeply rising background from decays B^ 1 —t J/i/jK ± at higher masses. Neither the B$ signal 
nor background components are constrained to positive values in order to obtain unbiased results for rates close to zero. 



VII. SIGNAL YIELD 

The measured decay rate of a particle resonance as 
a function of the invariant mass of the final state, 
is described by the relativistic Brcit-Wigncr function 
(RBW) [1] convoluted with detector resolution. 

To obtain the detector resolution, we use simulated 
— > J /i/jf '2(1525) decays where the J/tp is forced to 
decay into two muons and the /2(1525) into two kaons. 



The fitted M(K + K~) distribution for the simulation is 
shown in Fig. 15. Fixing the mass and natural width 
parameters at their input values, M = 1525 MeV and 
Tq = 73 MeV, and using the range parameter [1] R ~ 
5.0 GeV" 1 , we obtain a(M) = 22 ± 1 MeV. 

We fit the signal yield versus M(K + K~) from data, 
as shown in Fig. 10, to an incoherent sum of the J = 2 
component and a constant continuum term. The result 
is shown in Fig. 16. The fit yields 629±157 J^ 1525 ) 
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FIG. 9: Invariant mass distribution of the meson pair from 
the B® candidates in the mass range 5.4 < M(J /ipK + K~) < 
5.6 GeV, under the K Tt^ hypothesis. The curve shows the fit 
of a relativistic Breit-Wigner J = 2 resonance at a fixed mass 
of 1.43 GeV and with a floating width, over a background 
described by a second-order polynomial function. 
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FIG. 10: The B ( > signal yield as a function of M(K + K~) 
obtained from the fits shown in Fig. 8. The outer and inner 
error bars correspond to the statistical uncertainties with and 
without systematic uncertainties added in quadrature. 



FIG. 12: Angular dependence of the detection and event se- 
lection acceptance of the decay B° — > J /ipf 2(1525) from sim- 
ulations as a function of (a) cos Oh, (b) 4>h, and (c) cos ip. 
The acceptance is found to be independent of the angle <f>n- 
For cos Oh and cost/', we fit the acceptance dependence with 
symmetric fourth-order polynomial functions. 




FIG. 11: Definition of the decay angles Oh, 4>h, and tp in the 
helicity basis for the sequential decay — > J/ipX, J/ip — > 
/j+fj,-, X -> K+K-. 



events and 345±76 events for the constant term in the 
mass range 1.4 < M(K + K~) < 1.7 GeV. The fraction 
of the nonresonant term, assumed to be the S wave, in 
this mass range is 0.35 ± 0.09. 



VIII. RATIO OF B(B° -> J/i/>/2(1525)) AND 
B(B° -> J/ip<t>) 

To determine an absolute branching fraction for the 
B® — > J/^i/2(1525) decay, efficiencies, branching frac- 
tions, and the cross section need to be known, as well 
as the integrated luminosity. However, terms common 
to the B° -> J/^/2(1525) and the B° -> J/ipcj) branch- 
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as shown. The fits assume a two-Gaussian B® signal tem- 
plate and a background composed of a reflection of the B° — > 
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falling background from multibody decays of B mesons at 
lower masses and a steeply rising background from decays 
B ± -> J/ipK* at higher masses. These backgrounds are par- 
ticularly acute for cos(^i) > 0.8 making it impossible to obtain 
a reliable fit in this region. 



600- 



o 
c 

Uj' 



400- 



: D0 Run II, 10.4 fb" 1 / / 




+ Data / / 




— J=2 + J=0 / / 




J=2 / / 




J=1 / / 




j=o J/ 


■ . . . i . 





0.2 0.4 0.6 0.8 1 

|COS\|/| 

FIG. 14: The | cos ip\ distribution for the decay B° -> J/ipX, 
X -» K + R- in the range 1.45 < M(K + K~) < 1.6 GeV. The 
curves are best fits assuming pure J = (dashed line), pure 
J = 1 (dashed-dotted line), and pure J — 2 (solid line). 



> 

CD 



600- 



?40Q 

tn 

| 200 
lli 



■ D0, MC j 


- Simulation 


A 




/ 4 


- RBW, J=2 


1.4 1.6 


1.8 



M(K + K") (GeV) 

FIG. 15: The M(K + K~) distribution from the simulation of 
the decay B® — > J/ipf 2(1525) fitted by the relativistic Breit- 
Wigner function [1]. 



>400h 



: D0 Run II, 1 0.4 f b" 1 _i_ Data 

— Full Fit 




RBW, J=2 
Const 



1.4 1.6 1.8 2 

M(K + K") (GeV) 

FIG. 16: Fit to the B° yield versus M(K + K~) as obtained in 
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signal described by a relativistic Breit-Wigner and nonreso- 
nant constant term for the S wave. The outer and inner error 
bars correspond to the statistical uncertainties with and with- 
out systematic uncertainties added in quadrature. 



ing fractions cancel in their ratio. A measurement of 
the relative branching fraction Rf requires the yields 
of the two decays, ATb°->j/^(1525) and N B o^j/^, and 
the reconstruction efficiencies of the two decay modes, 

B^JfM , B°->J/^(1525). 
crcco dim crcco 



R 



fU4> 



B(B° S -> .7/^(1525); #(1525) -> K+K~) 



B(B® -> J/ip(f>; <t> -> K+K- ) 



N, 



B°->-J/i/>/2(1525) x £ reco 



N R o_ 



B°->J/^(1525) 



(4) 



where -/Vb°->.j/i/)/^(i525) is determined in the mass range 
1.4 < M(K + K~) < 1.7 GeV and N B °^j/^ in the mass 
range 1.01 < M(K + K~) < 1.03 GeV. The only differ- 
ence in the event selection for the two channels is the 
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M(Kn) > 1 GcV condition applied for the J/ty/^ (1525) 
candidates. 

The yield of the B® — > J/tp4> decay is determined by 
fitting the data, shown in Fig. 17, with a double Gaus- 
sian function for the signal, a second-order polynomial for 
background, and the reflection of the decay B® — > J/ipK* 
taken from simulations. The total number of B® — > J/ip4> 
events is 4064 ± 105. 




FIG. 17: Invariant mass distribution of B s candidates 
with ct > 200 fim for events in the mass range 1.01 < 
M(K + K~) < 1.03 GeV. A fit to a sum of a double Gaussian 
B° — > J/ip<j> signal (dashed-dotted line), a quadratic combi- 
natorial background (dotted line), and the reflection of the 
decay B° — > J/ipK* (892) (dashed line) is used to extract the 
B° yield. 



We use simulated samples of the two decay processes 
to determine the reconstruction efficiencies. For the de- 
cay B® — > J/ipj 2(1525) the efficiency is measured to 
be (0.122 ± 0.002)% and for the decay B° s -> J/^ it is 
(0.149 ±0.002)% (where the uncertainties are due to MC 
statistics), yielding Rf^/4, = 0.19 ± 0.05 (stat). 

The denominator in Eq. 4 may include a contribution 
from the K + K~ S wave, and no correction is made, al- 
lowing the ratio to be recalculated for different 5-wave 
fraction inputs. 

IX. SYSTEMATIC UNCERTAINTIES 

The main contributions to systematic uncertainties are 
summarized in Table I. They are evaluated as follows: 

• Kg (1430) width: Wc vary the Kq width within its 
uncertainty of 0.08 GeV [1]. 

• ^(1430) and if|(1430) templates: We vary the 
shape of the K$ and templates by altering the 
widths of the dominant Gaussian component within 
statistical uncertainties. 

• Combinatorial background shape: As an alterna- 
tive, we use a second-order polynomial to describe 



the combinatorial background. Wc also vary the fit- 
ting mass range from 5.25 — 5.70 GeV to 5.2 — 5.8 
GeV. 

• Signal shape: We vary the B® mass scale within its 
uncertainty in data of 0.08% and the width of the 
core Gaussian component by ±10%. 

• Trigger efficiency: Due to the mass difference be- 
tween the /2(1525) and <j> resonances, there is a 
small difference between average muon momenta in 
the two channels. Approximately 3% more J/ipfj) 
events have a leading muon with pt > 15 GeV 
and about 3% more 7/ ip<p events have both muons 
with pt > 3 GeV. We therefore estimate that there 
is approximately a 3% difference in the fraction of 
events that can be accepted by the trigger between 
the J/ifj(j) and 7/^/2(1520) signals. Trigger sim- 
ulations confirm this estimate. We apply the 3% 
correction to Rf and assign an absolute 3% sys- 
tematic uncertainty. 

• M{K + K~) dependence of efficiency: The 
M(K + K~ ) dependence of the efficiency for recon- 
structing the ^(1520) resonance is obtained from 
a simulation. Wc assign a 2% uncertainty due to 
the statistical precision of the MC sample. 

• Helicity dependence of efficiency: The B° s — > J/ip4> 
signal acceptance is obtained from a MC sample 
generated under the assumption that the final state 
is not polarized, i.e., with the final state distributed 
uniformly in helicity angle cos 6h . We compare this 
signal acceptance with distributions corresponding 
to pure helicity and 1 and assign a systematic 
uncertainty equal to the difference. 

• /2(1525) mass and natural width: The uncertainty 
on the mass of the /2 (1525) resonance of 5 McV [1] 
leads to the uncertainty on Rf^/4, of 3%, while the 
uncertainty on the natural width of 6 MeV leads to 
the uncertainty on Rfy<p of 0.7%. 

• Bg — > J/ip4> signal shape: The B® — > J/t/xj) signal 
yield is sensitive to the signal mass model because 
of the presence of the B® — > J/ipK*(890) reflec- 
tion that peaks near the signal. We assign the sys- 
tematic uncertainty as one half of the difference 
between fit results for single and double Gaussian 
distributions for the signal mass model. 



X. SUMMARY AND DISCUSSION 

We confirm the observation of the decay B® —> 
7/^/2(1525) previously observed by the LHCb Collab- 
oration [3] and measure the ratio of branching fractions 
of the decays B° -> 7/^/^(1525) and B° s ->■ 7/^0 to be 
R fi/rj) = 0.19 ± 0.05 (stat) ± 0.04 (syst). The fit to the 
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TABLE I: Sources of systematic relative uncertainty on R^y 



Source 


Uncertainty (%) 


#£(1430) width 


5 


#5(1430) and # 2 *(1430) templates 


10 


Combinatorial background shape 


10 


Signal shape 


12 


Trigger efficiency 


3 


M(K + K~) dependence of efficiency 


2 


Helicity dependence of efficiency 


3 


/a (1525) mass 


3 


/2(1525) natural width 


1 


-B3 — > J /4>4> signal shape 


4 


Total 


20 



background-subtracted signal — >• J/ipK + K , assum- 
ing an incoherent sum of the J = 2 resonance /2(1525) 



and a constant continuum term, assigns the fraction of 
0.35 ± 0.09 to the constant term. The fit to the helicity 
angle i/j in the K + K~ rest frame finds the K + K~ reso- 
nance to be consistent with J = 2 with a fit probability of 
0.27 and preferred over J = or J = 1, for which the fit 
probabilities are 2.8 x 10" 2 and 9.8 x 10 -3 , respectively. 
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